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ABSTRACT 

The  Apollo  Lunar  Program  utilized  efficient,  i.e.,  Earth-return,  transearth 
trajectories  which  employed  parking  orbits  in  order  to  minimize  energy  requirements. 
This  thesis  concentrates  on  a  different  type  of  transearth  trajectory.  These  are  direct- 
ascent,  hyperbolic  trajectories  which  omit  the  parking  orbits  in  order  to  achieve  short 
flight  times  to  and  from  a  future  lunar  base.  The  object  of  the  thesis  is  the  development 
of  a  three-dimensional  transearth  trajectory  model  and  associated  computer  program  for 
exploring  trade-offs  between  night-lime  and  energy,  given  various  mission  constraints. 
The  program  also  targets  the  Moon  with  a  hyperbolic  trajectory,  which  can  with  a  time- 
reversed  trajectory;  be  used  for  targeting  Earth  impact  points.  The  first-order  model  is 
based  on  an  Earth-centered  conic  and  a  massless  spherical  Moon,  using  MathCAD 
version  3.0.  This  model  is  intended  as  the  basis  for  future  patched-conic  formulation  for 
the  design  of  fast  Earth-return  trajectories.  Applications  include  placing  nuclear- 
deterrent  arsenals  on  the  Moon,  various  space  support  related  activities  and  finally 
protection  against  Earth-threatening  asteroids  and  comets  using  lunar  bases. 
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I.      INTRODUCTION 

This  chapter  provides  a  general  description  of  this  thesis  and  is  divided  into  two 

sections.    The  first  section  presents  background  information  pertaining  to  translunar 

trajectory  design.    The  second  section  delineates  the  objectives  and  limitations  of  this 
project. 

A.   BACKGROUND 

While  the  impending  mission  of  establishing  a  permanent  lunar  base  stands  in  our 
future,  it  is  appropriate  to  study  this  type  of  mission  by  investigating  various  translunar 
trajectories.  The  trade-off  between  flight  time  and  energy  required  must  be  determined 
prior  to  considering  the  type  of  translunar  trajectory  to  design.  The  energy  required  for 
the  hyperbolic  trajectory  can  then  be  compared  to  that  of  the  Apollo  mission  to 
determine  the  increase  in  energy  and  thereby  cost  increases. 

The  Apollo  Lunar  Program  of  the  1960s  utilized  energy-efficient  translunar 
trajectories.  This  thesis  will  analyze  a  different  class  of  trajectories  named  fast  Earth- 
return  trajectories.  These  trajectories  utilize  a  direct-ascent  from  the  lunar  surface  and 
omit  the  lunar  parking  orbit  used  in  the  Apollo  missions.  These  direct-ascent  trajectories 
differ  considerably  from  translunar  trajectories  previously  considered  and  discussed  in  this 
thesis.  (Wadsworth,  1991) 


B.       OBJECTIVES  AND  LIMITATIONS 

The  objective  of  this  thesis  is  to  develop  a  trajectory  model  and  associated 
computer  program  for  systems  studies  of  fast  Earth-return  trajectories  from  a  lunar  base. 
The  scope  of  the  project  involves  a  step  process  which  evolves  in  complexity,  as  the 
design  moves  from  a  massless  point-Moon,  to  a  massless  spherical  Moon  intercept 
solution.  Finally,  the  groundwork  is  laid  for  a  patched-conic  model  which  culminates  in 
a  solution  for  intercept  on  a  Moon  with  finite  lunar  mass.  A  novel  feature  is  the  use  of 
a  time-reversed  trajectory  for  targeting  Earth  impact  points.  In  other  words,  the  lunar 
impact  point  serves  as  the  lunar  launch  point  for  a  time-reversed  trajectory.  The 
complete  patched  conic  approximation,  as  discussed  in  Chapter  II,  is  not  developed  in  the 
program  due  to  scope  of  work  and  time  constraints. 


II.      ANALYSIS  APPROACH 


This  chapter  sets  the  foundation  of  the  thesis  hy  delineating  a  flow  chart  of  tasks 
that  are  required  to  be  accomplished  to  reach  the  objective  of  the  thesis.  A  flow  chart 
represented  in  Figure  1,  summarizes  the  work  to  be  accomplished.  The  last  block  of  the 
diagram  represents  further  research  work. 


DESIGN  A  LUNAR  TRAJECTORY 
TO  INTERCEPT  A  MASSLESS 
POINT-MOON  MOVING  IN  A 
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PERTURB  THE  INITIAL  CONDITIONS  AT  THE  EARTH  LAUNCH  POINT  BY 

ADDING  A  MANEUVER  AT  THE  SPHERE-OF-INFLUENCE  CROSSING  OF 
THE  EARTH  AND  MOON 
THE  COMPOSITE  PATCHED  CONIC  WILL  USE  A  TRAJECTORY 
THAT  CAN  BE  ADJUSTED  TO  INTERSECT  THE  SPHERICAL 

MOON  AT  ANY  DESIRED  SURFACE  IMPACT  POINT  ON  THE 
NEAR-SIDE  HEMISPHERE  OF  THE  MOON. 


Figure  1.   Translunar  Trajectory  Design  Flowchart. 

Precision  translunar/transearth  trajectory  design  requires  time-consuming,  step-by- 
step  numerical  solution  of  the  exact  differential  equations  of  motion,  including  lunar  and 
solar  ephemerides.     For  preliminary  design,  the  computation  time  can  be  reduced  a 


hundred-fold  or  more  by  employing  a  closed-form  "patched-conic"  approximation  to  the 
precision  trajectory.  Since  the  patched-conic  is  a  composite  of  the  two  conies,  it  can  be 
formulated  explicitly  in  terms  of  elementary  functions  of  the  key  parameters.  This 
provides  better  insight  into  design  trade-offs  than  numerical  computation.  The  patched 
conic  approximates  the  gravitational  effects  of  the  Earth-Moon  system.  It  ignores  the 
small  perturbation  due  to  the  Sun. 

Two  approaches  were  considered  for  developing  the  patched-conic  model.  Both 
approaches  employ  Earth-centered  and  Moon-centered  coordinate  frames  with  respect 
to  which  the  conic  trajectories  are  formulated.  Although  the  conic  trajectories  are  planar 
with  respect  to  the  Earth-centered  and  Moon-centered  frames,  they  form  three- 
dimensional  curves  as  viewed  in  an  inertial  barycentric  frame.  In  the  first  approach,  the 
coordinate  frames  are  non-inertial,  being  attached  to  the  Earth  and  Moon  which  travel 
in  (idealized)  circular  orbits  (with  respect  to  the  Earth-Moon  inertial  frame).  The  second 
approach  employs  non-rotating,  inertial  Earth-centered  and  Moon-centered  coordinate 
frames.  These  frames  are  attached,  respectively,  to  a  fictitious  Earth  and  fictitious  Moon 
which  travel  at  constant  velocity  with  respect  to  the  barycentric  frame.  The  implications 
of  the  second  approach  are  briefly  stated  in  Appendix  A.  The  advantage  of  the  first 
approach,  pursued  in  this  thesis  research,  is  described  in  section  A.  (Wadsworth,  1991) 

A.       TRANSLUNAR  TRAJECTORY  DESIGN 

The  first  consideration  given  to  designing  a  translunar  trajectory  is  to  decide  on 
what  translunar  trajectory  types  will  be  used.  Since  this  thesis  concentrates  on  fast  Earth- 


return  trajectories,  eccentricities  greater  than  one  will  be  utilized.   The  reason  for  using 
non-inertial  coordinate  frames  is  discussed  in  section  A. 2. 

1.  Translunar  Trajectory  Types 

Figure  2  illustrates  all  possible  trajectory  types  that  can  be  considered  for 
Earth-Moon  trajectory  design.  The  radius  of  the  Earth  (R^),  radius  of  the  Moon  (RM), 
spacecraft  true  anomaly  at  launch  (fE),  and  spacecraft  true  anomaly  at  moon  intercept 
(fM)  are  projected  in  Figure  2. 

Since  only  trajectories  with  eccentricities  greater  than  1.0  are  being 
considered,  only  type  I  and  III  hyperbolic  trajectories  will  be  used  in  thesis  research  (see 
Table  1). 

2.  Non-inertial  Coordinate  Frames 

The  Earth-centered  and  Moon-centered  coordinate  frames  are  non-inertial 
frames,  since  they  are  in  a  state  of  gravitational  free-fall  about  the  Earth-Moon 
barycenter  (center  of  gravity).  As  a  consequence,  the  composite  patched-conic  provides 
a  better  approximation  to  the  restricted,  three-body  spacecraft  motion  than  would  be 
suggested  by  simply  joining  two  conies. 

Within  the  Moon's  region-of-influence,  the  terrestrial  perturbing  acceleration 
of  the  spacecraft  is  nearly  the  same  as  the  centripetal  acceleration  of  the  Moon-centered 
coordinate  frame.  In  this  frame,  these  accelerations  nearly  cancel;  this  can  be  deduced 
from  the  facts  that  the  distances  of  the  spacecraft  to  the  Earth  and  the  Moon  to  the 
barycenter  are  comparable  and  the  mass  of  the  Earth  is  comparable  to  the  equivalent 
mass  at  the  barycenter.    For  the  Earth-centered  phase,  there  is  a  similar  benefit  over 
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Figure   2.       Equatorial    plane   projection. 


TRANSLUNAR  TRAJECTORY  TYPES: 

TYPE  1  posigrade  :  arc  ab 

TYPE  I  retrograde     arc  dc 

TYPE  II  posigrade  :  arc  abc  (  e  <  1  ) 

TYPE  II  retrograde  :  arc  deb  (  e  <  1  ) 

TYPE  III  posigrade  :  arc  dab 

TYPE  III  retrograde  :  arc  adc 

TYPE  IV  posigrade  :  arc  dabc  (  e  <  1  ) 

TYPE  IV  retograde  :  arc  adeb  (  e  <  1  ) 


Table  I.   Trajectory  Types. 


most  of  the  trajectory,  especially  for  parabolic  and  hyperbolic  trajectories.  This 
cancellation  does  not  occur  in  the  case  of  patched-conics  based  on  inertial  frames  as 
described  in  Appendix  A  (Wadsworth,  1965). 

The  Earth-centered  translunar  conic  trajectory  may  be  either  elliptic  or 
hyperbolic.  The  latter  corresponds  to  shorter  flight  times  and  is  the  only  case  that  needs 
consideration  for  fast  trajectories.  For  the  same  reason,  the  moon-centered  conic 
approach  trajectory  is  also  modeled  as  hyperbolic  in  this  study.  Both  the  Earth-centered 
and  Moon-centered  frames  are  non-rotating.  The  motions  of  the  Earth  and  Moon  about 
their  common  barycenter  are  approximated  by  coplanar  circular  orbits.  In  general,  the 
plane  of  these  circular  orbits  does  not  coincide  with  the  three  dimensional  trajectory  of 
the  spacecraft.  (Wadsworth,  1991) 

B.       SOLUTION  FOR  POINT-MOON  INTERCEPT 

The  least  complicated  form  of  a  first  order  solution  is  for  a  point-moon  intercept. 
Figure  3  illustrates  the  relevant  parameters  of  the  problem.  Appendix  B,  utilizing 
MathCAD  3.0,  numerically  solves  the  first  order  solution  of  point-Moon  and  massless 
spherical  Moon  intercept.    In  Figure  3  the  following  assumptions  are  made: 


•  The  Earth-centered  cartesian  coordinate  frame  is  non-rotating  (indicated  by 
X,Y,and  Z). 

•  The  Moon-centered  cartesian  coordinate  frame  is  non-rotating  (indicated  by  x,y,and 
z  with  the  x-axis  parallel  to  the  X-axis,  and  the  z-axis  tilted  at  angle  iM  for  the  Z- 
axis.) 


The  initial  point-Moon  intercept  solution  represents  a  simple  three-dimensional 
model;  three-dimensional  considerations  are  incorporated  in  the  first-order  solution  solved 
in  Appendix  B.   These  considerations  will  he  discussed  in  section  B.2. 

1.  Three  Dimensional  Point-Moon  Model 

Figure  3  illustrates  the  three-dimensional  model  of  the  first-order  solution  at 
point-moon  intercept  while  Figure  4  further  illustrates  aspects  of  the  three-dimensional 
model.  Figure  5  delineates  a  flowchart  of  tasks  to  be  accomplished  in  the  three- 
dimensional  model  of  the  first-order  solution  of  point-Moon  and  massless  spherical  Moon 
intercepts.  This  flowchart  summarizes  steps  used  to  construct  Appendix  B.  Symbols  and 
abbreviations  are  defined  in  the  Tabic  of  Symbols  and  Abbreviations.  The  hyperbolic 
trajectory  yielded  an  energy  requirement  of  63.2  km2/sec2  while  the  Apollo  mission 
required  a  minimum  energy  requirement  of  61.8  kirr/sec2.  This  results  in  a  2%  increase 
in  energy  by  using  the  hyperbolic  trajectory,  and  as  a  consequence,  will  increase  launch 
costs.  Considering  that  the  time  of  flight  is  reduced  by  more  than  50%,  the  2%  energy 
increase  penalty  seems  feasible  for  this  mission. 

2.  Three  Dimensional  Considerations 

Figure  6  illustrates  the  relationship  of  orbit  eccentricity  and  spacecraft  true 
anomaly  at  moon  intercept.  Figure  7  is  an  exploded  view  of  Figure  6  for  user 
interpretation  which  shows  eccentricities  between  1.0  and  2.0.  These  diagrams  must  be 
referenced  by  the  user  to  insure  appropriate  inputs  for  solution  in  Appendix  B.  Figure 
8  illustrates  a  cone  sweeping  out  a  circular  area  of  the  Earth.  The  Earth  surface  within 
the  cone  defines  a  zone  of  inaccessibility  for  Earth  launch  sites,  given  the  specified  true 
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Figure  5.   Translunar  Trajectory  Program  Flowchart. 

anomaly  at  moon  intercept.  This  fact  constrains  inputs  to  the  program  in  Appendix  B. 
If  the  launch  site  is  found  to  be  in  the  restricted  cone  region  a  different  trajectory  must 
be  used. 
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Figure  8.   Three  dimensional  considerations. 

C.       SOLUTION  FOR  MASSLESS  SPHERICAL  MOON  INTERCEPT 

The  final  solution  in  Appendix  B,  is  the  massless  spherical  Moon  intercept.  An 
iterative  technique  is  employed  to  obtain  an  impact  on  the  lunar  surface.  This  technique 
is  illustrated  and  described  in  section  C.l.  Selenographic  coordinates  are  utilized  to 
obtain  the  solution  for  the  massless  spherical  moon  intercept  and  are  described  in  section 
C.2. 

1.        Iterative  Technique  for  Lunar  Surface  Impact 

Figure  9  illustrates  the  method  to  ascertain  the  lunar  surface  impact  point 
iteratively.   The  iterative  process  is  defined  in  the  following  steps: 
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•  Determine  point-moon  intercept  coordinates  (R  =  pEM,  fM,  $M,  and  tM  are  inputs). 

•  Select  trial  value  of  radial  distance  :  R0  «  pEM-D. 

•  Calculate  the  corresponding  flight  time  from  perigee,  t0,  and  separation  flight  time 
D0. 

•  Find  the  next  iteration  :  R,  =  R0  +  (D„  -  D)cosA.0. 

•  Rn«Dn-rRnDn  =  cos A n and  Rn  =  R„ .,  +  (D„  ,-D)cosA.n.,  (bold  print  represents  vector 
quantities). 

•  Continue  until  n  =  10  or  |  DIV,-D  |  <  1  km. 

This  iterative  technique  is  employed  in  Appendix  B,  however,  solutions  are  calculated 
manually  by  the  program  user  without  a  terminal  loop. 
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2.        Selenographic  Coordinates 

Impact  points  on  the  near-side  of  the  moon  are  expressed  in  selenographic 
coordinates,  for  which  the  zero  degree  longitude  and  latitude  points  at  the  Earth.  The 
North,  South,  East,  and  West  directions  indicated  on  Figure  13  for  latitude  and  longitude, 
are  for  an  observer  on  the  Earth.  The  lunar  librations  (Roy,  1988)  are  neglected  in  this 
definition. 

a.  Euler  angles 

Figures  10  and  11  express  Euler  angle  transformations  which  are 
incorporated  in  Appendix  B  to  relate  selenographic  coordinates  to  Earth-centered 
coordinate  systems.  The  first  Euler  angle, ((f)),  is  zero  because  this  angle  is  not  required 
in  this  analysis.  The  second  Euler  angle  rotation, (0),  is  the  angle  iM  in  Appendix  B.  The 
third  Euler  angle  rotation,(TM),  is  the  angle  $M  in  Appendix  B.  (Goldstein,  1951) 

b.  Lunar  Impact 

Figure  12  illustrates  a  selenographic  coordinate  solution  for  lunar  impact 
of  the  spacecraft.  Since  time  reversal  is  employed  in  Appendix  B,  the  Moon  moves 
backwards  in  its  orbit  for  the  translunar  trajectory.  The  sign  on  wM  is  reversed  in 
Appendix  B  for  this  reason. 

The  selenographic  latitude  for  surface  impact  is  designated  8 
and  the  selenographic  longitude  for  surface  impact  is  designated  a.  Figure  13  illustrates 
these  latitudes,  longitudes  and  their  signs.  The  solution  for  a  has  a  factor  of  ir  added  to 
it  to  insure  that  the  positive  X-axis  points  to  the  Earth.    Figure  14  displays  a  locus  of 
impact  points  for  varying  time  offsets  (toff).   The  axes  units  are  in  radians.    Varying  toff 
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is  equivalent  to  varying  the  input  value  of  the  lunar  sweep  angle  (4>M).  By  providing  an 
offset  from  the  point-Moon  intercept  declination  (<5M  ofl),  the  user  can  target  from  North 
to  South  (see  Figure  13).  This  factor  can  he  comhined  with  (toff)  to  produce  an  arbitrary 
location  of  the  impact  point  on  the  lunar  near  side. 
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Figure    14.       Locus    of    lunar    impact   points. 

D.       PATCHED  CONIC  APPROXIMATION 

The  method  of  calculation  for  the  massless  spherical  moon  intercept  model  can  also 
be  used  for  sphere-of-influence  calculations  if  the  radius  of  the  sphere-of-influence  is 
substituted  for  the  radius  of  the  Moon.    The  patched  conic  lunar  impact  point,  using  a 
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finite  lunar  mass  model,  can  then  be  approximated  by  an  iteration  process  identical  to  the 
process  used  in  Appendix  B  (Bate  et  al.,  1971). 

From  the  final  patched  conic  iteration,  the  impact  point  for  the  Moon-centered 
conic  trajectory  can  be  calculated.  This  calculation  is  made  using  the  same  method  as 
Appendix  B. 
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[II.      APPLICATIONS 

Fast  Earth-return  trajectories  have  numerous  applications,  such  as  for  requirements 
in  evacuating  a  lunar  base  in  an  emergency  or  the  establishment  of  a  lunar-based 
strategic  deterrent  force.  This  chapter  will  discuss  such  applications,  their  development, 
and  technical  feasibility. 

A.       STRATEGIC  PEACE  INITIATIVE  CONCEPT 

The  ABM  treaty  prohibits  nuclear  weapons  in  space.  If  we  look  at  a  future  world, 
this  treaty  might  possibly  be  altered  or  not  exist.  The  Strategic  Peace  Initiative  is  a 
hypothetical  plan  to  place  nuclear  weapons  on  the  Moon  and  to  eliminate  them  from  the 
Earth.  The  reasoning  for  this  replacement  is  to  prolong  the  strike  capability  and  reaction 
time  that  now  exists,  thereby  increasing  stability. 

1.        Background  Information 

This  section  describes  a  hypothetical  Strategic  Peace  Initiative,  a  plan  to 
eliminate  the  risk  of  global  nuclear  war  and  simultaneously  promote  peaceful  uses  of 
space.  The  SPI  concept  could  be  a  step  toward  reducing  present  nuclear  confrontation 
as  well  as  "denuclearizing"  the  Earth.  (Wadsworth.  1988) 

Rather  than  using  our  current  posture  of  SDI,  the  SPI  concept  can  encourage 
peaceful  exploitation  of  space  (Gray,  1985).  SPI  resources  could  provide  the  logistics 
support  for  a  new  era  of  nonmilitary  space  activities,  such  as  manned  space  exploration, 
servicing  geostationary  communication  satellites,  safe  disposal  of  radioactive  waste,  clean- 
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up  of  hazardous  space-debris,  and  eventually,  deflection  of  Earth-threatening  asteroids 
or  comets. 

Although  the  basing  of  deterrent  strike  forces  in  space,  particularly  on  the 
Moon,  has  already  been  proposed;  this  proposal  is  quite  unique.  The  novelty  of  an  SPI 
concept  is  that  by  a  new  treaty,  the  United  States  and  Soviet  Union  would  share  the 
"high  ground"  of  space.  As  with  both  the  INF  and  Start  treaties,  success  of  the  SPI  treaty 
would  depend  on  whether  it  serves  the  mutual  interests  of  parties  involved.  (Wadsworth, 
1988) 

The  free  world's  offensive  deterrent  policies  have  prevented  a  global  nuclear 
confrontation  to  date.  Unfortunately,  this  policy  is  dangerously  unstable  because  of  the 
critically  short  decision  time  for  launching  a  counterstrike  to  a  perceived  nuclear  attack. 
SPI  would  eliminate  or  greatly  reduce  the  danger  of  accidental  nuclear  war  by  increasing 
decision  time  from  minutes,  in  the  case  o(  an  SLBM  attack,  to  a  more  comfortable  two 
days. 

The  primary  justification  for  SPI  is  survival  in  a  nuclear  age,  consequently  the 
emphasis  in  this  proposal  is  on  deterrence.  A  secondary,  economically-compelling 
justification  is  that  much  of  the  SPI  investment  in  global  security  could  also  serve  the 
civilian  space  program.   They  both  require  similar  logistics  elements: 

•  A  heavy-lift  launch  vehicle 

•  A  space  tug  (space/lunar  transfer  vehicle) 

•  Earth/Moon  space  communications  and  tracking  networks 

•  A  manned  lunar  base 
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2.        SPI  Architecture 

The  minimum-energy  transearth  strike  trajectory  has  a  4.7  day  flight  time. 
Doubling  the  injection  energy  yields  a  three  day  trajectory.  High-energy  strike  trajectories, 
taking  two  days  or  less,  still  permit  nuclear  strike  weapons  to  escape  the  Earth-Moon 
system  if  aborted.  Appendix  B  provides  such  a  time  of  flight  (approximately  1.9  days). 
Shortening  the  flight  time  from  4.7  days  to  2  days  results  in  only  a  two  percent  increase 
in  transearth  trajectory  energy,  which  is  equivalent  to  providing  the  additional  velocity 
increment  typical  of  an  ICBM.  Atmospheric  re-entry  velocity  does  not  increase 
significantly  between  4.7  day  and  2  day  trajectories,  which  simplifies  re-entry  vehicle 
design.  To  minimize  vulnerability  to  intercept,  the  trajectories  must  be  of  the  direct- 
ascent  type  without  lunar  parking  orbits,  in  contrast  to  the  more  economical  supply 
shuttle  trajectories  which  must  also  be  utilized.  (Wadsworth,  1988) 

The  SPI  communications  networks  and  surveillance  systems  would  preferably 
be  located  at  the  natural  synodic  satellite  points,  designated  LI  through  L5,  which  are 
stationary  with  respect  to  the  Earth-Moon  system  and  would  minimize  station-keeping 
fuel.  These  points  are  therefore  an  ideal  location  for  placing  these  communications 
satellite  systems.  (Farquar.  1970) 

The  two  opposing  lunar  base  zones,  each  approximately  one  hundred 
kilometers  in  diameter  would  contain  the  strike  force  missile  silos  and  manned 
command/support  complexes.  These  zones  would  be  located  on  approximately  opposite 
sides  of  the  moon.  They  might  occupy  polar  sites  to  take  advantage  of  full-time  sunlight 
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for  solar  power.    Determining  optimum  zone  locations  requires  analysis  if  the  optimum 
strike  trajectory  missions  for  both  bases.  (Wadsworth,  1988) 

Similar  to  the  deployment  of  the  stages  of  SDI,  SPI  would  be  implemented 
in  stages  and  be  subject  to  treaty  compliance  verification  (possibly  a  United  Nations  role) 
at  each  stage.    Both  superpowers  would  deploy  the  following  five  stages: 


•  NON-NUCLEAR  DEPLOYMENT  :  Lunar  bases,  communications,  and 
surveillance  systems  deployed. 

•  FLIGHT  TESTS  :  Test  strike  force  with  dummy  warheads  to  demonstrate  accuracy 
and  reliability. 

•  NUCLEAR  DEPLOYMENT  :  Equip  lunar  bases  with  nuclear  missiles  and  phase 
out  Earth-based  strategic  nuclear  forces. 

•  VERIFICATION  :  Surveillance  of  lunar  bases  and  Earth-Moon  space;  inspection 
and  monitoring  of  Earth-based  nuclear  device  manufacturing  and  test  facilities;  and 
inspection  of  translunar  cargos. 

.  STRATEGIC  REDUCTION  EVOLUTION  :  Bilateral  reduction  of  strategic 
nuclear  forces  to  minimum  required  for  asteroid  deflections;  increase  in  cooperative 
nonmilitary  space  ventures  by  utilization  of  lunar  bases. 


B.       SPI  TECHNICAL  FEASIBILITY 

Initial  deployment  of  SPI  relies  on  technology  that  has  or  will  soon  be  proven  by 
other  programs.  In  fact,  SPI  can  virtually  be  used  as  a  spin-off  of  developed  SDI 
technology.  Programs  such  as  TDRSS,  NASA  space  station,  OMV,  NASP,  HEDI,  and 
the  APOLLO  lunar  programs  are  examples  of  some  of  these  programs. 

Challenges  in  technology  i.\o  exist.  These  challenges  include  a  precision  guided  re- 
entry vehicle,  developing  "smart"  surveillance  techniques,  and  prevention  of  adverse 
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physiological  effects  of  long-term  residency  in  low  gravity.  At  an  interval  of  deployment 
of  ten  to  twenty  years,  the  fully  deployed  SPI  deterrent  force  system  probably  would  cost 
comparably  to  the  proposed  SDI  system.  (Wadsworth,  1988) 

Under  the  SPI  treaty,  a  sovereign  lunar  base  zone  would  be  assigned  to  each 
superpower  for  basing  its  nuclear  strike  force.  These  zones  would  be  located  on 
approximately  opposite  sides  of  the  moon.  The  optimum  size  and  location  of  the  zones 
would  be  a  compromise  among  the  following  six  basic  requirements: 


•  Minimize  Surveillance  Cost:  Requires  minimizing  zone  size;  locate  bases  on  the 
Moon's  visible  side  to  take  advantage  of  surveillance  redundancy  provided  by  Earth- 
based  sensors  to  back  up  the  spaced-bascd  sensors. 

•  Maximize  Nuclear  Survivability:  Requires  dispersal  of  hardened  missile  silos  within 
a  circular  area  of  at  least  one  hundred  kilometers  diameter. 

•  Minimize  Base  Vulnerability:  Requires  maximizing  separation  between  bases  to 
deter  pre-emptive  strikes  and  reduce  vulnerability  of  transearth  strike  forces  to 
interception. 

•  Minimize  Transportation  Costs:  Requires  selection  of  base  zone  locations  to 
minimize  fuel  cost  for  both  supply  shuttles  which  utilize  parking  orbits  and  strike 
forces  which  require  direct-ascent  trajectories.  (This  issue  will  be  difficult  to  resolve, 
because  if  both  strike  forces  were  to  take  advantage  of  the  Moon's  orbital  velocity, 
then  locations  near  the  lunar  poles  will  be  favored.  This  induces  a  strategic 
problem,  because  the  energy  requirements  for  striking  Northern-hemisphere  Earth 
targets  will  differ  between  the  polar  bases.). 

•  Avoid  Base  Overflight:  Requires  designing  supply  shuttle  parking  orbits  not  to 
overtly  opposing  base. 

•  Minimize  Communication  Costs  :  Requires  designing  base  communications  on  the 
near-side  of  the  Moon.  Far-side  Moon  bases  require  more  communication  system 
elements  because  its  orientation  with  the  Earth  does  not  meet  full-time  line-of-sight 
communications  link  requirements. 
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C.  SPACE  SUPPORT  ACTIVITIES 

A  large  variety  of  space  support  programs  can  be  served  by  SPI.  For  example, 
lunar-based  astronomy  would  be  free  of  restriction  imposed  by  Earth's  atmosphere.  A 
lunar  base  could  be  used  as  a  staging  platform  for  a  manned  expedition  to  Mars,  taking 
advantage  of  a  smaller  round  trip  propulsion  requirement  for  liftoff  from  the  Moon, 
rather  than  from  the  Earth.  The  servicing  of  geostationary  satellites  could  be 
accomplished  with  SPI  resources.  Transport  of  concentrated  nuclear  waste  to  a  safe 
dumping  site  in  the  Moon  would  be  unhampered  by  groundwater  problems  as  on  the 
Earth  (Rosen,  1981).  SPI  resources  could  be  used  for  debris  removal  from  low-Earth 
orbits.  This  is  an  increasingly  serious  hazard  to  space  missions.  Also,  safe  disposal  of 
space  junk  or  dead  satellites  in  higher  orbits  could  be  an  SPI  priority. 

Using  the  SPI  surveillance  network,  the  nuclear  arsenal  of  the  Moon  can  be  used 
to  deflect  Earth  threatening  asteroids  or  comets.  Collision  with  a  ten-kilometer  diameter 
asteroid  has  the  explosive  force  of  ten  million  one-megaton  nuclear  bombs.  With  an 
explosion  of  such  magnitude,  the  impact  dust  cloud  will  plunge  the  world  into  a  deep 
freeze  similar  to  "nuclear  winter"  but  much  more  devastating.  Fossil  records  indicate  that 
such  an  asteroid  could  account  for  extinction  of  the  dinosaurs  and  cause  for  the  Ice  Age 
(Malove,  1985). 

D.  SUMMARY 

The  Strategic  Peace  Initiative  Concept  is  a  revolutionary  plan  for  strategic  nuclear 
deterrence.    It  is  a  plan  which  provides  double  value;  a  solution  to  accidental  global 
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nuclear  war  by  increasing  reaction  time  lo  at  least  two  days,  and  promotion  of  peaceful 
uses  of  space  by  developing  permanently  manned  lunar  bases  and  their  infrastructures. 
The  uniqueness  of  SPI  is  that  it  requires  bilateral  treaty  agreement  between  the 
United  States  and  the  Soviet  Union.  Although  the  politics  of  such  a  treaty  requires 
extensive  thought,  verification  of  the  treaty  can  be  accomplished.  By  looking  to  the 
future  of  peace,  perhaps  we  can  be  at  peace  with  ourselves  and  with  future  generations. 
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IV.      CONCLUSIONS  AND  RECOMMENDATIONS 

This  final  chapter  summarizes  the  findings  if  this  thesis.  It  also  recommends 
research  for  further  work  concerning  this  topic. 

A.  FINDINGS 

This  thesis  utilizes  the  selection  of  an  efficient  and  operationally  meaningful  set  of 
inputs  and  terminal  parameters,  to  lead  to  a  transearth  or  translunar  trajectory  targeting 
solution,  which  rapidly  converges  in  an  iteration  process. 

B.  RECOMMENDATIONS 

This  thesis  lays  the  groundwork  for  an  accurate  patched  conic  method  of  designing 
fast  transearth  trajectories.  Recommendations  include  the  finalization  of  the  patched 
conic  approximation  and  the  iterated  impact  point  for  the  Moon-centered  conic 
trajectory.  A  development  of  the  SPI  concept  to  include  a  cost/mission  analysis  would 
also  be  very  useful. 
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APPENDIX  A.      PATCIIED-CONIC  BASED  ON  INERTIAL  FRAMES 

The  patched  conic  trajectory  can  be  formulated  in  terms  of  Earth-centered  and 
Moon-centered  coordinate  frames  which  arc  inertial,  being  attached  to  a  fictitious  Earth 
and  Moon  which  travel  at  constant  velocity  with  respect  to  the  barycentric  frame.  The 
Earth-centered  conic  trajectory  approximation  implies  the  lunar  gravitational  force  is 
ignored  so  that  the  Earth  is  not  accelerating  with  respect  to  the  barycentric  frame. 
Consequently,  after  the  instant  the  spacecraft  is  launched,  the  center  of  gravitational 
force,  the  Earth,  is  replaced  by  a  fictitious  Earth  which  is  imagined  as  moving  in  a  straight 
line  at  constant  velocity  in  the  barycentric  frame.  In  the  meantime,  the  massless  Moon 
continues  on  its  circular  orbit  in  the  barycentric  frame.  In  terms  of  the  translating,  Earth- 
centered  frame,  the  lunar  orbit  appears  as  a  compound  curve,  rather  than  a  circle.  This 
fact  must  be  taken  into  account  in  targeting  the  conic  trajectory  for  lunar  intercept  or 
crossing  of  the  lunar-sphere-of-influence.  The  same  considerations  apply  to  the  Moon- 
centered  conic  phase  where  the  Moon  is  approximated  by  a  fictitious  Moon  moving  in  a 
straight  line  at  constant  velocity  in  the  barycentric  frame.  For  greatest  accuracy,  the  line 
is  taken  tangent  to  the  Moon's  circular  orbit  in  the  barycentric  frame  at  the  time  of 
intercept.   This  requires  an  iterative  solution.  (Wadsworth,  1991) 
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APPENDIX  B.      LUNAR  TRAJECTORY  PROGRAM 


INPUT  PARAMETKltb1. 


S/C  declination  at  earth  launcb  point 


(input  in  degrees) 


5t. 


E 


180 


■28.3 


6E  =  0.493928178314395 

radiens 


Irae  anomaly  of  s/c  el  moon  intercept 
(input  in  degrees) 


M 


164 


180 
fM  =  2.8623399732707 

radians 


inclination  of  tin:  moon  orbit  plane 
(input  in  degrees: 

20.5  >  i  >  10.!)) 


lM  • 


180 


20 


iM  =  0.349065850398866 

radians 


eccentricity  of  tfce  IranH&Gar  trajectcrj 

(eocer  >  =  eccen     ) 


mm 


geocentric  altitude  of  launch  point 
(altitude  in  kilometers) 


eccen  :     1.01 


hv   -  1 
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lunar  weep  angle  al  moon  intercept 
(input  in  degrees) 


'M 


180 


•170 


0M  =  2.96705972839036 

radioes 


posigrede  or  retrograde  orbit  determination 
poagrade  =  +-1 
retrograde  =  -1 


aE  OBT    -  1 


5M    off    :=fl 


not  iterated  in  program 


IN1ERKEDMTE  OUTPUT: 


SM  :=  asin[sin[iM]'sinj>M]] 

5M  =0.059426145347819 
(spacecraft  deeiiflalioii  at  moon  iolercept) 


Ml 


cos   <f>M 


_cos[oMj_ 


aM1   =-0.986549223316463 


:M2  ': 


[cos[iM]-sin[^M] 


cos 


^rvr 


«M2  =0.163464460888851 


-  ,.  f     f[1_0!Ml] 
aj4  ■-  2*atan J 

aM2 


aM  =  2.97739131857372 
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INPUT  CONSTANTS: 

pem  =  384400  km  :  eaKh-nooc  meae  distance 

rE  :=  6378  km  :  earth  meai  radwis 

rM  :=  1738  km  :  moon  mean  radios 

u  E  :=  7.292-10-5    rad/sec :  earth  angchr rotation  rale 


0)m  ■=  1.525*10     •— 
M  180 


cj  M  =  2.661627109291352*10         rad/ttC  : 

moon  neon  orbilal  rotational  rale 

2     3 
k.E  :=  398600  km  /sec  :  earth  gravitational  constant 

2     3 
kM  :=  4903  km  /sec  :  moon  gravitational  cot^toot 

OUTPUT  YARUBIES  and  FORMULAS: 
w  M  :=  -i  'w  m  (die  to  Ike  root  moving  backwards  in  its 

orbit  for  Uc  IransJoear  trajectory) 

RE1  :=  rE*hE 

R  ei  =  6.379*10         km  :  (s/c  radial  distance  at  earth  lainch) 
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RM1     "  PEM 


RM1  =  3.844* io5  km  :  (s/c  radial  distance  al  mooi  inlercepl) 


r  := (radial  distance  ralio} 

RE1 


r  =  60.26022887599937 


fp  :=  acos 


r*[  1  +ccccn*cos[  f  j^j  "1  "J  -  1 


eccen 


fc  =  0.726288833059664 


fEM  -  fM"fE 


fEM  =  2.136051140211037 


s'gn fEM  '=  sin[fEM] 


NOTE:  for  S 


fem 


signf£M  =  0.84445214187218 


if  sign  is  positive  Ibee  use  +1 
if  sign  is  negative  use  -1 


Sfcm    -  * 


_[C0S[{EM]"sin[gE]'sin[5M]1 
ttEM1  [co.[8E]-cos[«M]] 


a  EM  1   =  "0-64 J  4527643620 1 7 
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aEM2    -J1_aEMl    *sfem'aE  OBT 


aEM2  =0.767162532383019 
[1_aEMl] 


"EM       2'alan 


!EM2 


NOTE: 
»f  <*EM1         =    "I   ll»ei> 

IS  ir"«E_OBT 

if  «EM1     =    +'    ^eB 

is  0 


a 


EM 


EM 


aEM  =  2.267186781944015 


«E    -  «M"«EM 


a£  =  0.710204536629705     ( s/c  rigbl  ascecsioi  al  earth  launch  ) 


5pi    -  _x  +  f \A  ~d 


El 


M"°M 


6  El  =  -0^38678825666912 


5E2   "  *_fM~5M 


8E2  =  0.219826534971273 


fM_MAXl    =  s'n[5Ei,sin[5M]  +  cos[aEM],cos[5E]'cos[5M 


fM_MAXl  =-0-53563101113312 


M_MAX2 


1-f 


M_MAX1 


^M_MAX2=0-84445214I87218 
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M_MAX 


:=  2*atan 


[1-fM_MAXl 


M_MAX2 


fM_MAX  =  2.136051140211037 


If: 
Ihee 


6E1  <6E<5E2 

fM<fM_MAX 


is  constrained 
otherwise  use 


«E 

as  calcololed 


Aa  :=  asin 


tan[*E]'tan[*M]- 


cos 


[[Ml 


[cos[5E]-cos[5M] 


Aa  =  1.570796326794897  -  0.496321 2320603 15i 


Cecal  rami: 


aEMIN    =  aM"2_Aa 


aEMIN  =-0-164201335016073  +  0.496321232060315i 


aE      m«st  be  greater  Uac 


aE 


M1N 


ccccn 


mm 


1- 

1 

r 

1- 

2 
r 

"cos[fEM]  + 

1 

~2 

LN 

r 

ccccn^j,  =0.974647175723671 
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cncrgymm 


_'[-kE*[l+ccccnmin]" 


[2'RMl] 


( =  minimim  specific  energy 
al  b«Dob  ) 


energy  min  =  -1.023795999276086 


#E1 


sin[3M]-sin[5E]-cos[fEM] 
cos[8Ej-sin[fEMJ 


j8Ei  =  0.421410460300705 


0E2 


cos[5M]-Sin[aEM] 

sid[[em] 


/3E2  =  0.906870014913464 


0E  :=  2*atan 


0E2 


digit  path  azinolt)  at  leoicb: 


*>/S>-x 


0E  =  1.135796261808675 


cos[5El*sinr/SE  1 


jj  =0.798478510707822 


«2  :=  V1"' 


i2  =0.602023311789352 


34 


i   -  2*atan 


1-i 


l-> 


(IraRstonar  trajectory  indinalioi: 
0<  =  i<=x) 


i  =  0.646032654504068 


aE 


0E1 


Ql 


aEQ!  =  0.699990269559789 


tan 


aEQ2   " 


[«*] 


Un(i) 


aEQ2  =  0.71415238O463452 


oEq  :=  2*atan 


1-aE 


Ql 


aE 


Q2 


uEq  ■  0.795412455434794 


Q  :=  ag-aEfj 


(Irasslunar  Irajeclory  longRvdo 
of-ascendiiiE-iiode) 


O  =  -0.085207918805089 


Cfrj  :=  cosf  Sjr  l'cosfaEQ  1 


CEI  =0.616325580024169 


35 


E2 


sin[«E]" 


C^  =  0.787491447196649 


Cg  :=  2*atan 


1-C 


El 


E2 


CE  =  0.906728172342672 


(sweep  angle  al  laoocfe: 

ar>c  >  =  o) 

E 


Rp  = 


Rjri  *T  1  +ccccn*cosrf£l 


1  +  eccen 


(perigee  distance) 


Rp  =  5^7010652841956*10 


V„  := 


IkE'(l  +cccen) 
R^ 


(perigee  velocity) 


Vp  =  11.99319583652338 


E    : 


2-k, 


KE1 


(l-CCC«3l)'k.g 


R, 


(s/c  velocity  at  lac b eh) 


VF  =  11 2\  107362995628 
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R: 


Rp*(i    ccccn)        (apogee  distance;  valid  only  for  eccea  (  1 


1  -  ccccn 


ignore  for  eccec  >  =  l) 


Ra  =  5.57010652841956*10 


^asy       walan 


u 


ecccn   _ 1 


fasy  =  3.000756780023376 


(hyp^rtoBc  asymptotic  line  anomaly 
balf-angle :  valid  oily  for  eccen^  1 
ignore  for  eccen  <  -  1} 


D=     2-t^y-K 


(deflection  asgie) 


D  =  2.859920906456958 


■  asy  ' ~      p 


,.p 


ecccn 


(directrix- to- locos  distance  or  axial  ofsel 

of  asymptotic  Irom  locos ;  valid  only  for 
eccen  >  1  igiore  for  cceerj  <  =  1) 


<*asy  =  1.108506348725081  MO 


VE1 


VVP 
RE1*VE 


YEl  =  0.934111286367402 


YE2 


eccen'y^j  'sin^f  £ 1 
1  +  ccccn  •  cosf  f  £  1 


yE2  =  0356981938872315 


37 


yp  '=  2*alan 


1-yEi 

YE2 


(fligfet  pelh  angle  a(  IgidcIi; 


IT  -IT 

2       L     2 


yE  =  0365034946775507 


2*kE 


M 


R 


Ml 


( 1  _ccccn)*k.£ 
R„ 


(s/c  velocity  al  moot  intercept] 


VM  =  1.670175762283258 


?M1 


RP*VP 

*mi'vm 


yM1  =0.104052584194718 


fMl    =cos[fM] 


lM2  :=sin[fM] 


?M2  := 


eccen-yM1  -f^ 
l  +  [ecccn-fM1] 


YM2  "  0.99457179716821 


yM    =  2-atan 


>M2 


(fligkl  paid  angle  al  moot  iilercept) 


yM  =  1.466555060009538 


38 


0M1  :='2'[cos[aM-Q]] 


0M1  =-0.600145987662734 


0M2  '- 


cos 


[«*] 


0M2  =0.799890488437212 


Pm  •=  2'aUn 


0M2 


(fligkl  patfc  azimcU  al  noon  iolercepl) 


0M  =  2.214479932657055 


XM  :=  RMl  'co»[*m]"«»[«M  ]  (RCCTANGUUR  COORDINATES  OP 

THE  MOON  AT  POINT-MOON 


XM  =  -3.785601002578928 '10 


INTCRCEPT) 


YM  :=  RMl,cos[5M]"s'n[«M] 


YM  =  6.2724820252415%  MO 


Zm  -  Rmi  ••"[*!«] 


ZM  =2.282996752157727*10 


E     E  (Iroislcoar  Irajeclory  ai*g«meit- 

u  =  0.180439339283008  of- perigee) 

39 


CM  :=  fM+w 


Cj^  =  3.042779312553709        (sweep  angle  al  moon  intercept) 


■M* 


M 


[[rM2-sin[5M]-i2-cos[CM-yM]]-cos[aM-Q]] 


cos[aM] 

+  -i*[ii*r'Mrsin[aM~Q]] 


XM.  =  -1.665737958076395 


M* 


M 


cos 


i  •[[yM2-*in[3M]*'2,cos[cM-rM]],sin[aM-Q]] 

;['i'YMrcos[aM-Q]] 


YM-  =  0.048513058379955 


zw  -Vm^'^Cm-ym] 


ZM-  =  -0.00545766973131 


[Xm"%]+[ym-Ym'Hzm,zm'] 

RM* ~ 

RM1 


(radial  *e)ocRy  at  po'ml-moon  mterce^l) 


RM-  =  1.648023680688986 


40 


2**'R 


P  .5 


L5 

(circular  oitll  period  for  radios  R  ) 

P 


Tp  =  4.137199589247244-10 


1  "COS 


El 


M 


"smjTjTj 


fE1  =037999703209726 


l-cos[fM] 

IT"' 


Ml 


sm^tM] 


iff- 


fM1  =  7.115369722384208 


Ehl  :=  In 


eccen  -  1  , 

1+  *fEl 

^  1  +  ccccn 


cccen  - 1  _ 
1  +  eccen 


Ehl  =  0.05361869041727 


Eh2  :=  ^ 


cccen  -  1  , 
1  + eccen 


cccen  -  1  . 
/\|  1  +  eccen 


Ml 


En2  =  1.103630544045282 


41 


shl 


-rl 


eccen  -  1  , 

'  h *fEi 

1 +  eccen 


rl 


1  + 


ccccn  - 1 
1  +  eccen 


El 


Eshl  =0.053644386077804 


sh2 


eccen-  1 


>Jl+eccen     M1 


-rl 


leccen  -  1   , 

1+  ' *fMl 

1+ eccen     M1 


Esh2  =  134171391602585 


f cccca  *Eshl]-Ehl 

*E1  := [T 

2"ir,(eccen  -  1) 


tE1  =0.089467283524084 


feccen'EB^l-E^ 

tM1  :=~ — ; ^TT 


2*ir*(eccen-l) 


tM1  =40.02754953820076 


tE:ssTp'*El 


(earth  lacncl  Digkl  lime  relative 
to  perigee  passage) 


tE  =  370.144008646907 
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tM^Tp^Ml 


(poiil-moon  intercept  fligtl  line 

relative  lo  perigee  passBge) 


1M  =  1.656019615080179-10 


*EM  -  ^"'E 


(poiit-moon  Inaidsnar  Irajedory 
flight  lime) 


tEM  =  1.65231817499371 '10 


INPUTS: 


R0  :=  382986^976625   (aiy  radial  distance) 


OUTPUTS: 


Vn  := 


2-kE]   [(l-eccen)*kE] 


Ri 


R, 


Vn  =  1.6724659380983 


_[Rp'(l+eccen)] 
~R^ 


R,  =  0.029233197611653 


In  '■-  2'atan 


[  ccccn  +  1  -  R  x  1 
f  ecccn  -  1  +  Rx  1 


f0  =  2.861954051456415 


43 


Cn    -  w  +fr 


C  0  =  3.042393390739423 


5q  :-  asin[sin(i)*sinrCo]^ 


50  =  0.05%5775 1937721 


Aa 


Ml 


[K8m]] 

[sin(i)-tan[CM]] 


^Ml  =-0.996881642139323 


Aa 


M2 


_tan[5M] 


tan(i) 


AaM2  =0.078911289214964 


Aa^  :=  2'atan 


[*-*»Ml] 


Aa 


M2 


Aa  M  =  3.062599237378809 


tan 


Aaj 


[*o] 


[sin(i)-tan[C0]] 


Aa!  =-0.996857191524672 
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tan 


A£K2 


[M 


tan(i) 


Aa2  =  0.079219566431173 


Aot0  :=  2*atan 


[1-Aa,] 


Aa 


6a  n  =  3.062289992046826 


ai  :=  aM_A«M 


a t  =  -0.085207918805089 


a0  :=  Aa0  +  at 


a0  =  2.977082073241737 


T,-[Vvp] 


yj  =  0.104293660037493 


[ccccn'r^sinjfo]" 
f  1  +ccccn*cos[  f0] 


y2  =  0.994546546158593 
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■Yq  '■-  2*atan 


I'-.r 


>2 


y0  =  1.466312665340157 
Pi  =  -0.600131267822715 


P  2- 


cos 


[»o]J 


02  =  0.799901532303508 


0O  :=  2*atan 


*2 


0O  =  2.214461530464985 


X0   -  Ro*cos[50],cos[ao] 


X0  =  -3.771433240671579*10 


Y0   -  Ro'cos[60]'sin[a0] 


(RECTANCUIAR  CO0RDJNATCS  OF 
THE  SPACECRAFT  ON  ITS  TRAJECTORY) 


Y0  =  6.260990806241751  *10 


Z0   =  Ro'sin[50] 


Z0  =  2.283455104046146  MO 
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_[l-cos[f0]] 

HfoJ 


MO 


sin 


fM0  =  7.105421102050813 


h3  :=ln 

1  + 

|(ccccn  -  1) 
\(  1  + ccccn) 

*fM0 

1- 

|(  ccccn  -  1 ) 
aJ  ( 1  +  ccccn ) 

*fM0 

Ey  =  1.101755458701495 


sh3 


1 


( ccccn  -  1 ) 
( 1  +  cccen ) 


MO 


1  + 


( cccen  -  1 ) 
( 1  +  ccccn ) 


1-1 


*fM0 


Esh3  =  1J38578545032142 


MO 


[[ccccn'Esh3]-Eh3] 
[2"ir-(ccccn-l)1"5] 


tM0  =  39.8219787493874 


*0    -  Tp"!MO 


trt  =  1.647514741 24978*  10 
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tnff  =0  (lime  oVsel  in^al) 


Atq  =  »M_to+ton 


Ar0  =  850.4873830398719 


♦mi  -  cos[«m]'cos[8m] 


^M1  =-0.984807753012208 


sin 


♦  M2  :" 


[5m] 


K'm] 


sin 


♦  M2  =  0.17364817766693 


^M  :=  2*atan 


['-♦mi] 

♦M2 


♦  M  =  2.%705972839036 


♦  MO  -♦M"a)M*ArO 


♦  MO  =  2.%9323408665169 
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XM0  =  pem*cos[*mo]  (RECTANGULAR  COORDINATES  OF THE 


XM0  =-3.787 10231 6833942  MO5 


YMO  :=PEM'cos['M]"sin[*M0] 


YM0  =  6. 19 1940094679928  *104 


ZM0  :=PEM*sin['M]"sin[*M0] 


ZM0  =  2.253681886822944  *104 


NOON  AT  mCl  ) 

0 


49 


DM0  :=  4[X0-XM0]2+[Y0-YM0]2+[Z0-ZM0]2] 


DISTANCE  BETWEEN  THE  SPACECRAFT  AND  THE  LONAR  CENTER  AT  TIME  t 

0 


DM0  =  1.738000000264939  *103  km 


'  THE  COAL  IS  TO  DRNE  THIS  DISTANCE  TO  EQUAL  THE  RADIUS 
Of  THE  MOON  FOR  AN  IMPACT  ON  Tffi  SURFACE  OF  THE  MOON 


error  -  DMO"rM 


error  =  2.649394446052611  -10  7 


A       .    [[XmO-Xo]-Xo+[Ymo-Yo]-Yq+[ZMo-Zo]-Z0] 


ApQ  =  0.81263979792074  '  cos[Xn]     ' 


Rncxt    -  R0  +  DMO"ApO 


Rnext  =  3.843986656315015-105  (  Mnber  l°  <*  *>"**  "^^  °f 

R0) 
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X  -  X\m  - X 


M    AM0 


X  =  150.1314235014295 


Y = Yu-Y 


M"  YM0 


Y  =  805.4193056166769 


Z-ZxmZ 


M~^M0 


Z  =  293.148653347831 


INPUTS: 

0  (EHERANGICNPUTS) 


$  :=  i 


M 


*M    =  *M0 


*M  =  2.969323408665169 
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OUTPUTS: 

(THE  MATRIX  STRUCTURE  IS  NUMBERED  AS  IUlOWS:| 


A] 

Bl     C," 

A2 

D2    C2 

A3 

b3  c3 

(  ROTATION  M  TRANSLATION  ) 


A1  :=  cos[*M]*a>s(*)-cos(0)*sin(*)'sin[*M] 


A  j  -  -0.985198313432347 


Bt  :=  -sin[*Ml*co8(*)-cos(0)*sin(*)*cos[¥Ml 


B  i  =  -0.171418444777836 


Cj  :=  sin(0)'sin(^) 


Cj  =  0 


A2  :=  cos[*M]*sin(*)*cos(0)'cos(*)»sin[*M] 


A->  =  0.161080647624348 


B2   =  -sinr*Ml*sin(^)+cos(0)*cos(*)*cos[*Ml 


B  ->  =  -O.925783585143099 


C2   =  -sin(0)*cos($) 


C2  =  -0342020143325669 
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A3  :=  sin(<?)*sin[*M] 
A3  =  0.058628561051585 
B3  :=  sin(0)*cos[*M  1 

B  3  =  -0336957668364338 


C3  :=  cos(0) 


C  %  =  0.939692620785908 


A    = 


Al 

Bj    Cj 

A2 

B2    C2 

A3 

B3    C3 

B 


C  :=  A'B 


C  = 


-285.9729519600736 

-821.7234494824876 

12.87940451383525 
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C0  =  -285.9729519600736 


Cj  =  -821.7234494824876 


C2  =  12.87940451383525 


xlun   ~  c0 


Ylun    -  cl 


zlun  :_  c2 


d  '■-  alan 


lun 


4Xlun   +Ylun  J 


i  =  0.014801756974441 


SEIENDCRAPBJC  LATITUDE  FOR  SURFACE 
IMPACT  ON  THE  KOON 


a  •-  -acos 


lun 


4XI"n   +Y»un  j 


+  * 


a  =  1235890035212915 


SEIEWCRAPBIC  LQNCITUDE  FDR  SURFACE 
IMPACT  ON  THE  MOON 
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